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Semiconductor quantum dots coherently driven by pulsed laser are fundamental physical systems 
which allow studying the dynamical properties of confined quantum states. These systems are 
attractive candidates for a solid-state qubit, which open the possibility for several investigations in 
quantum information processing. In this work we study the effects of a specific decoherence process, 
the spontaneous emission of excitonic states, in a quantum dot molecule. We model our system 
considering a three-level Hamiltonian and solve the corresponding master equation in the Lindblad 
form. Our results show that the spontaneous emission associated with the direct exciton helps 
to build up a robust indirect exciton state. This robustness against decoherence allows potential 
applications in quantum memories and quantum gate architectures. We further investigate several 
regimes of physical parameters, showing that this process is easily controlled by tuning of external 
fields. 

PACS numbers: 73.21.La, 73.40.Gk, 03.65.Yz 



I. INTRODUCTION 

The advance on the manipulation and dynamical con- 
trol of quantum states under the action of coherent radi- 
ation has recently become a subject of intense research in 
condensed matter physics. Dynamical control is a neces- 
sary step for the implementation of any protocol associ- 
ated with Quantum Information Processing (QIP) [l]. In 
this sense, semiconductor quantum dots (QDs) has been 
proved to be an ideal candidate. Using strong resonant 
laser pulses and different probe techniques, several differ- 
ent groups have successfully demonstrated coherent ma- 
nipulation of the exciton population of a sing le QD. [H3 
They demonstrated a process known as Rabi oscillation 
which is indeed a proof of the exciton qubit rotation. Un- 
like atoms, however, QDs suffer from unavoidable varia- 
tion in their size, and the presence of a surrounding envi- 
ronment with which they may interact strongly, making 
the entire system to lose its phase quickly |8Hig|. The 
main interest in QDs arises from their characteristic dis- 
crete energy spectrum, and its great flexibility in change 
it, not only by manipulation of their geometric structure, 
but also with the application of external gates. A natural 
next step for the development of such system is to put two 
quantum dots together and allowing them to couple. A 
lot of work has been done in this direction, where beauti- 
ful examples of a molecule formation have been achieved 
pl] - |22j . However, the coherent dynamics of such objects 
under strong laser pumping remains largely unexplored 
experimentally, and our work can give further insight to 
help the experimental development. 

In this paper, we study the effects of the spontaneous 
decay in the excitonic states of a self-assembled semicon- 
ductor quantum dot molecule (QDM) coupled by tun- 
neling and under the influence of an external electromag- 
netic field. We use a standard density matrix approach in 
the Lindblad form to describe the system dynamics and 
our results indicate that the spontaneous decay of the 
direct exciton helps to build up a coherent population of 
the indirect exciton (electron and hole in different dots). 



which has a longer lifetime due to its spatial separation 
with small overlap of the wave function. This effect is 
robust is robust to the changes of external parameters 
and in order to describe it we describe the physical sys- 
tem and the detailed theoretical model in Sec. [TTl then 
in Sec. Illll we show the results of numerical calculations, 
followed by our conclusions in Sec. IIVI 



II. THEORY AND MODEL 

The physical system we consider here is an asymmet- 
ric double quantum dot coupled by tunneling. Electrons 
and holes can be confined in cither dot and we can use a 
near-resonant optical pulse to promote electrons from the 
valence to the conduction band, creating an electron-hole 
correlated state known as exciton. Electrons or holes can 
then tunnel to the other dot, creating an indirect exciton. 
An external electric field, applied in the growth direction, 
brings the individual levels of electrons or holes into reso- 
nance, favoring the tunneling. Nevertheless, in asymmet- 
ric QDM structures it is even possible to control which 
type of carrier, electron or holes, tunnels ^23i] . In this 
situation, we can safely neglect the tunneling of holes as 
the electric field brings one level (conduction band) more 
into resonance while makes the other (valence band) more 
out of resonance. With this assumption, the dynamics of 
the QDM can be modeled by a simple three-level sys- 
tem, where the ground state |0) is a molecule without 
any excitation, |1) is the system with one exciton in the 
left dot, while |2) is the system with one indirect exciton, 
after the electron has tunneled. The schematic configu- 
ration of levels and physical parameters are shown in the 
Fig. [1] where we also include the decoherence channels 
associated with spontaneous emission of excitonic states 

(rj, r§). 

Using the rotating wave approximation and dipole ap- 
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FIG. 1: Scheme of energy levels with physical parameters on 
Hamiltonian 



proximation, the system Hamiltonian is written as 24l |: 



m = E;=o^i-ij>oi +T^'=(ii)(2i + i2)(ii) 

+;il7(e'"^*|0)(l| + e-'"^*|l)(0|), 



(1) 



where ujj are the frequencies of |j)-th states {j = 0, 1, 2), 
Te is the tunneling coupling, is the frequency of the 
applied laser and the dipole coupling is 51 — ■ 

l^\l)/2h, where ~jt is the electric dipole moment and ~^ 
is the amplitude of incident field. The intensity of in- 
cident field can be easily controlled to provide available 
conditions for the coherent control of the system quan- 
tum state. We also assume a low-intensity incident pulse, 
so the Rabi frequency is significatively smaller than the 
intraband excitation energy [1^, $7 ^ wio ~ uJi ~ ujq 
and uj ~ wio- Under this assumption, we might consider 
that only the ground-state exciton can be formed in our 
system. 

Applying the unitary transformation 



U — exp 



(|1)(1|-|0)(0| + |2)(2|) 



(2) 



and using the Baker- Hausdorff lemma [2f| , we obtain a 
time-independent version of Hamiltonian ^ written as 
follows: 



, , -5i 2hn 
2 ' 2Te 62 



(3) 



where 5i = hiwio — ojl) is the detuning between the 
frequency of optical pulse and exciton transition, 82 = 
(5i + 2fiW2i and Uij is the optical transition between i and 
j energy states. 

To taking into account the effects of decoherence, 
we used the Liouville-Von Neumman-Lindblad equation 
given by [23|: 



dp{t) 

dt 



■[H,p{t)]+L{p{t)). 



(4) 



Here, p (t) is the density matrix operator. The Liouville 
operator, L{p), describes the dissipative process. Assum- 
ing the Markovian approximation, Liouville operator can 
be written as psl j: 

m ^lY.^)i2\j)m){j\~m{^\ - \^)m, (5) 



where F* corresponds to the decoherence rates due spon- 
taneous decay from the state \i) to the state \ In order 
to investigate the dynamics associated with this physi- 
cal system, we solve the master equation and found 
the density matrix coefficients at certain time t. Writing 
Eq.Q in the basis defined by |0), |1) and |2) states, we 
obtain a set of nine coupled linear differential equations 
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In order to solve the set of equations © , we rewrite as 
p — Ap, considering p as a column vector and A being 
a square matrix associated with the coefficients of the 
coupled system above. The solution can be written as 



(7) 



where, Xj and Sij are the eigenvalues and the matrix 
formed by the eigenvectors of matrix A, respectively. 
Pij{0) are the elements of the density matrix operator 
at i = 0. 



III. RESULTS AND DISCUSSION 

For our calculations, we consider the following val- 
ues of physical parameters: fiwio — 1.6 eV [1, [29j . 
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n 0.05 - 1.0 meVj^, 113, ~ 0.33-6.6 /LteV [Mlp, 
and Fq ~ 10~''rj!) [S^j- The tunneling coupling, which 
depends the barrier characteristics and the external elec- 
tric field, was selected as: Te ~ 0.01 — 0.1 meV '33'] or 
Te ~ 1 — 10 meV [sj], for weak and strong tunneling 
regime, respectively. The system dynamics depends also 
from the detunings Si and 62- Experimentally, Si is con- 
trolled by varying the frequency of external laser. The 
value of S2 is changed by varying Si and the frequency 
transition CJ21, which can be done by manipulation of 
external electric field that changes the effective confine- 
ment potential. By varying this set of parameters we 
are able to perform a coherent manipulation of the wave 
function of the system. For all simulations, we consider 
|4'(0)) = |0) as initial condition. 
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FIG. 2; Dynamics of populations. Pi, of level \i} (i = 0, 1, 2) as 
function of time for different choices of spontaneous emission 
rate FJ with parameters 5i = 0, CJ21 = 0, f2=50 /leV and 
Ts=10 /ieV. We use gray line for Po, light-gray line for Pi 
and black line for P2. (a) Non-dissipative dynamics (Fq — 0); 
(b)FS=1.0 ^eV; (c)FJ=2.0 fxeY; (d)Fj=4.4' ^teV. 

Our first task is to analyze the effect on population dy- 
namics of the decoherence process associated with spon- 
taneous emission of direct exciton. In Fig. [21 we plot 
the probability of occupation associated with each of the 
three levels considering different values of spontaneous 
emission rate Fj. For the physical parameters consid- 
ered here, the non-dissipative dynamics (Fj = 0) shows 
that there are Rabi oscillations between the three levels of 
the system. This is illustrated in Fig. [2] a). The popula- 
tion of indirect exciton, state |2), depends directly on the 
parameter Tg, although the value of the coupling fl and 
detunings 61 and UJ21, has important effects on dynam- 
ics [l^]. The situation changes when spontaneous emis- 
sion is taken into account. As we expected, the Rabi os- 
cillations become damped. This can be seen in Figs.[2]b), 
c) and d). For long times and values of Fj high enough. 



as shown in the inset of Fig. [2]d), the Rabi oscillations 
are suppressed and the electron wave function tends to 
an asymptotic state. 

Now we focuss our attention on the formation of a sta- 
tionary state with high population of indirect excitonic 
level, 1 2). With a lifetime significatively longer (about 
lO'' times the direct exciton) [32], this particular state 
shows more potential for quantum information process- 
ing than the direct exciton, |1) state. In order to study 
the effects of several physical parameters on Hamiltonian 
(HI) and the decoherence, we study the behavior of aver- 
age occupation of state |2), defined as 

P2 = — / P2 (t) dt. 

''oo Jo 
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FIG. 3: (Color online) Average population of state |2), P2, 
as function of detuning, 5i and frequency uu2i for 0=50 /leV, 
Te=W fieV and Fg f» 10~*rJ. a)rj=4.4 ^eV. Inset: K for 
Fq = and the same values of Q and Te. b)ro=0.66 /ieV. 

In Fig. |31 we plot our results for P2, as function of 
laser detuning Si and frequency a;2i, considering two dif- 
ferent values of the direct exciton spontaneous emission 
rate, Fq. Bright colors are associated with high values 
of P2 , which means an efficient transference of the elec- 
tron from the first to the second dot. From our results, 
it is possible to conclude that a large occupation prob- 
ability of 1 2) is obtained if the detuning 61 is balanced 
with the applied electric field so that -I- a;2i — 0. We 
will named this condition as balanced detuning. The be- 
havior considering full resonance between the three lev- 
els (^1 ~ W21 — 0) deserves more attention. Let us de- 
fine an area associated with the full resonance condition 
^i|il'^2i| ^ 50/ieV: when spontaneous emission is not 
considered (Fq=0), the average population P2 is near to 
zero, as shown in Ref. ^2^ and in the inset of Fig. [3] a). 
Thus, full resonance condition is not a good experimental 
choice for an optimal creation of indirect excitonic state. 
Considering the effects of spontaneous emission Fg, we 
can observe a different behavior: the values of P2 at point 
{uj2i,Si) = (0, 0) increase from 0.05 (for Fj^ = O) to ~ 0.6 
(for Fj=4.4 fieV) and ~ 0.2 (for Fj=0.66 ^leY). This 
shows that for realistic direct excitons, with a non-zero 
spontaneous emission rate, the transfer of the electron 
between dots is more efficient. 
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FIG. 4: (Color online). Average population of state |2), P2, 
as function of coupling parameter Q, and frequency 0121 for 
resonant and non-resonant condition, a) FJ = and Si — 0; 
b) Fj = and 5i=100 peV; c) FJ=0.66 /xeV and 5i = 0; d) 
Fj=0.66 /xeV and (5i=100 fieV; e) Fj=4.4 /xeV and 5i = 0; 
f) FJ=4.4 neV and (5i = 100 /^eV. In all cases, Te = 10 ^eV and 
Fg « 10-*FJ 



In Fig. m we show our results for average population, 
P2, as a function of both, frequency UI21 and dipole cou- 
pling r2, for different choices of Fq considering Si — 0, 
Figs.lHa), c) and e), and (5i=100 ^eV, Figs.|l]b), d) and 
f). For all cases, we are able to populate the indirect 
exciton state, evidenced by bright regions with values of 
P2 larger than 0.3. At resonance condition, Figs.|4]a), c) 
and e), this bright area have a V-like form, with higher 
values of P2 concentrated on a small area associated with 
low values of and 0^21. For non-resonant condition, the 
symmetry between negative and positive values of u)2i is 
broken. Still, the large values of P2 are obtained when 
the condition Si + W21 — is fulfilled. The action of 
spontaneous emission can be analyzed by comparing the 
different situations shown in Fig. U] Higher values of pa- 
rameter rj are connected with higher values of average 
population P2. For example, in Fig.|3]a) when Fq = the 
maximum value of P2 — 0.36. Considering decoherence, 
the maximum value of P2 goes from 0.6 for Fq=0.66 fieY 
[Fig. He)] to ~ 0.8 for Fj=4.4 fieV [Figg]d)]. Also, the 
total area for highly efficient population of |2) state in- 



crease as spontaneous emission increase: both, the arms 
of the characteristic V area and the region with best val- 
ues of P2 become progressively large when the value of 
Tq increase. 

It is useful to check the combined effect of both, the 
tunneling and decoherence. It is expected a good trans- 
fer of population associated with higher values of Tg pa- 
rameter. This can be verified by comparing the results 
P2 without the effect of decoherence process with Te=10 
/ieV, Fig.|3]a), with the results considering a higher value 
of tunneling parameter Tg ~ 50/icV, Fig.[5]a). The effect 
of decoherence process is illustrated by Fig. [S]b). Notice 
that area on FiglS] with high P2 increase by the action 
of decoherence and the maximum value of P2 goes from 
~ 0.4, in Fig. [Ha), to ~ 0.8 in Fig.[5]b). 
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FIG. 5: (Color online). Average occupation of state |2), P2, 
as function of coupling parameter, Q, and frequency aj2i con- 
sidering Te=50 fieY, (5i = and Fg « 10-*Fj. (a) Fj = 0. 
(b) Fj=4.4 /ieV. 

After our analysis of P2, it is necessary to check the 
actual behavior of level population P2 . In FiglHl we plot 
P2 considering some choices of physical parameters as- 
sociated with our previous analysis (Figs. [2] - [5]). In 
all cases, we limit ourselves to full resonance condition 
{5i ~ UJ21 — 0). When dynamics is associated with sta- 
tionary states, the value of P2 depends on two aspects: 
the final value of P2 at stationary state and the time 
needed to reach this maximum value. In Fig. [B] a) we 
plot P2 for ri=50 fieV and Te=10 /xeV considering dif- 
ferent values of Fj. We can conclude that a higher spon- 
taneous emission rate of the direct exciton is connected 
with a faster evolution to the asymptotic value of P2. 
That means, the broadening effects (short lifetime) on 
the direct exciton are advantageous if we are interested 
on manipulate electronic wave function in order to cre- 
ate an asymptotic state with high values of P2 at short 
times. 

Next, we verify that the exact maximum value of P2 
is related with coupling parameters and Tg and, also, 
with the balanced detuning condition Si +(^21 — 0. From 
our calculations of average occupation of indirect exci- 
tonic state, we analise the behavior of P2 considering a 
set of Q and Tg parameters associated with maximum 
values of P2 for two different Fj rates. 

The evolution of population P2 is shown in Fig. [5] b) , 
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FIG. 6; (Color online) Population of indirect exciton state, 
P2, as function of time at stationary regime. (a)P2 for dif- 
ferent values of spontaneous emission rate Fj considering full 
resonance condition (5i — 0, iU2i = 0) for coupling parame- 
ters Q—50 fieV and T^ — IO ^eV. (b) P2 associated with the 
physical parameters for the highest values of P2 founded in 
FiglHand FigH 



blue (red) lines represent the rj=4.4 /ieV (rQ=0.66 /ieV) 
situation. For a fixed value of F J , we can define the char- 
acteristic time to as the time at which the system reaches 
the asymptotic value of P2 (for example to ~ 14ns for Fj 
= 0.66/ieV). This characteristic time, which depends on 
the value of Fq, allows us to distinguish two dynami- 
cal regimes: 1) for long time, i 3> io the population is 
essentially independent of time and its maximum value 
depends directly on the fi/Tg rate value (the population 
increases when this rate increases). 2) at short times, 
t < to, the dynamics does not depend on the value of 
ri/Tg, being governed mainly by Fj. Comparing all cases 
plotted in Fig. [5] b) we can conclude that the condition 
to obtain an asymptotic state with large values of the 
occupation P2 , associated with short characteristic times 
to, is given by i?- ~ Thus, it is possible to obtain 

experimentally optimized values of P2, by adjusting ap- 
propriately the laser intensity £7 for fixed values of 
and rj, which in turn can be obtained through optical 



spectroscopy. 



IV. SUMMARY 

In this work, we use a standard density matrix ap- 
proach in the Lindblad form to model the dynamics of a 
Quantum Dot Molecule under the influences of externals 
electric and electromagnetic fields, and in the presence of 
spontaneous emission. By numerically solving the den- 
sity matrix we show that the spontaneous decay of the 
direct exciton helps to build up a coherent population 
of the indirect exciton, which should have important ap- 
plications in quantum information processing due to its 
longer coherence time. 

We further investigate the efficiency of creation of in- 
direct exciton state as function of physical parameters of 
our model. For weak spontaneous emission rate, the sys- 
tem presents a Rabi oscillation and in the opposite limit 
the system rapidly build up a stationary population of 
the indirect exciton. Our results shown that the popula- 
tion of the indirect exciton is strongly influenced by the 
spontaneous emission of the direct exciton. We demon- 
strate that the indirect exciton, which has a longer life- 
time, is robust against the spontaneous emission process. 
Finally, at maximum average population conditions, we 
determined a relation between the relevant parameters of 
the system which allows us to obtain large populations 
of indirect exciton P2 ~ 0.9. 
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